A graceful entrance to braneworld inflation 
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Positively-curved, oscillatory universes have recently been shown to have important consequences 
for the pre-inflationary dynamics of the early universe. In particular, they may allow a self- 
interacting scalar field to climb up its potential during a very large number of these cycles. The 
cycles are naturally broken when the potential reaches a critical value and the universe begins to 
inflate, thereby providing a 'graceful entrance' to early universe inflation. We study the dynamics 
of this behaviour within the context of braneworld scenarios which exhibit a bounce from a collaps- 
ing phase to an expanding one. The dynamics can be understood by studying a general class of 
braneworld models that are sourced by a scalar field with a constant potential. Within this context, 
we determine the conditions a given model must satisfy for a graceful entrance to be possible in 
principle. We consider the bouncing braneworld model proposed by Shtanov and Sahni and show 
that it exhibits the features needed to realise a graceful entrance to inflation for a wide region of 
parameter space. 



I. INTRODUCTION 

The inflationary scenario is presently the favoured 
model for large-scale structure formation in the universe 
El II II II U M- (For reviews, see 0, & H EjJ ) . Infla- 
tion arises whenever the universe undergoes a phase of 
accelerated expansion. Within the context of relativistic 
cosmology, a necessary and sufficient condition for accel- 
eration is a violation of the strong energy condition, i.e., 
p + 3p < 0. This is equivalent to the condition 7 < 2/3, 
where 7 = (p + p)/p defines the equation of state pa- 
rameter. In the simplest versions of the scenario, such 
a condition is realised by a scalar 'inflaton' field that is 
rolling sufficiently slowly down its self-interaction poten- 
tial, 2 < V @. 

Given the recent successes of inflation when confronted 
with observations of the Cosmic Microwave Background 
(CMB) radiation and high redshift surveys E3> a key 
question to address is how the initial conditions for suc- 
cessful inflation were established in the very early uni- 
verse. This question was recently examined within the 
context of Loop Quantum Cosmology (LQC), which is 
the application of Loop Quan tum Gravity to symmetric 
states. (For reviews, see E2, El [Tj] ) . In this scenario, 
semi-classical modifications to the standard Friedmann 
equations cause a minimally coupled scalar field to be- 
have as an effective phantom fluid below a critical value of 
the scale factor, a*, with an effective equation of state, 
7eff < E El El E El. This causes a collapsing, 
isotropic, spatially closed Friedmann-Robertson- Walker 
(FRW) universe to undergo a non-singular bounce into 
an expansionary phase. On the other hand, if the field 
satisfies the strong energy condition above a* (where clas- 
sical dynamics is recovered) , such a universe will eventu- 
ally recollapse. The overall result, therefore, is that the 
universe oscillates between cycles of contraction and ex- 
pansion E3. 

If the field's kinetic energy dominates its potential en- 
ergy during these cycles, it will effectively behave as a 
massless field with its value either increasing or decreas- 



ing monotonically with time. In principle, therefore, the 
field may gradually roll up its potential (assuming im- 
plicitly that it is evolving along a region of the potential 
that is increasing in magnitude). As a result, the poten- 
tial energy will gradually increase during each successive 
cycle, with the net result that it becomes progressively 
harder to satisfy the strong energy condition on scales 
above a*. Eventually, therefore, a cycle will be reached 
when the strong energy condition becomes violated dur- 
ing the expansionary phase of the cycle (a > a*), and this 
will initiate a phase of inflation. This implies that slow- 
roll inflation may be possible even if the field is initially 
located in a region of the potential that would not gen- 
erate accelerated expansion, such as a minimum. This 
effect has been termed a 'graceful entrance' to inflation 

E3- 

A key feature of the dynamics described above is that 
the universe is able to oscillate because the effective 
equation of state satisfies 7 e g < 2/3 for a < a* and 
7cfi > 2/3 for a > a*. Since these are rather weak 
conditions, it is important to investigate whether simi- 
lar effects are possible in other cosmological scenarios, 
such as the braneworld paradigm. This is the purpose 
of the present paper. The braneworld scenario is moti- 
vated by string/M-theory [2I E^] and has attracted 
considerable attention in recent years. (See [2! 0, E|| 
for reviews). Our observable four-dimensional universe 
is interpreted as a co-dimension one brane propagating 
in a five- (or higher-) dimensional 'bulk' space. A num- 
ber of boun cing; braneworld models have been developed 
to date IH El El IH HI E3. In the model proposed 
by Shtanov and Sahni (S-S), for example, the extra bulk 
dimension is timelike and this results in modifications to 
the effective four-dimensional Friedmann equation that 
induce a non-singular bounce |26( . 

The paper is organized as follows. In Section II, we 
outline how a phase space analysis for a field with a con- 
stant potential can yield valuable insight into the cos- 
mic dynamics that leads to a graceful entrance for infla- 
tion. We then proceed in Section III to determine the 
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necessary and sufficient conditions for graceful entrance 
when the Friedmann equation has an arbitrary depen- 
dence on the energy density of the universe. We consider 
the S-S braneworld scenario as an explicit example in 
Sections IV- VI and determine the regions of parameter 
space where a graceful entrance is possible. We conclude 
with a discussion in Section VII. 



II. THE PHASE SPACE DESCRIPTION OF 
GRACEFUL ENTRANCE TO INFLATION 

The oscillatory dynamics of isotropic universes sourced 
by a scalar field with a constant potential proves very 
useful when determining the necessary conditions for a 
graceful entrance to inflation for a given cosmological 
model. In general, the dynamics is determined by three 
evolution equations for the scale factor, the Hubble pa- 
rameter and the velocity of the scalar field, respectively. 
In addition, there is the Friedmann equation, which rep- 
resents a constraint. This implies that the dynamics can 
be expressed as a two-dimensional system of equations 
and presented on a two-dimensional phase space. 

In the LQC scenarios considered previously, it proved 
convenient to parametrize the phase space in terms of the 
scale factor and Hubble parameter The qualitative 
nature of the phase space that can result in a graceful en- 
trance to inflation is illustrated in Fig. ^ For a positive 
potential taking values in the range < V < V C rit, where 
V^ r it is some critical value that is determined by the pa- 
rameters of the model, the phase space contains a centre 
and a saddle equilibrium point. Two types of behavior 
are therefore possible. Trajectories that are sufficiently 
close to the centre encircle it and represent trajectories 
that undergo eternal oscillations. On the other hand, tra- 
jectories which pass above the saddle point represent ini- 
tially collapsing universes which evolve through a bounce 
into an eternal (de Sitter) inflationary era. This is illus- 
trated qualitatively in the top panel of Fig. ^ 

If the value of the potential is increased, the effect on 
the phase space is such that the centre and saddle point 
are positioned closer to one another. At the critical value 
V = V^rit, the points merge and disappear. The phase 
space for V > V CX it therefore contains no equilibrium 
points and all trajectories represent initially contract- 
ing universes which bounce and subsequently inflate, as 
shown in the lower panel of Fig. ^ 

This simplified phase space structure enables us to de- 
velop a graceful entrance mechanism for more realistic 
cosmologies, where the potential is field-dependent. Let 
us assume that the scalar field is initially located in a 
region of the potential such that V(4>) < V^rit and that 
the universe is undergoing oscillations about the centre 
equilibrium point. The dynamics in this regime will be 
dominated by the field's kinetic energy and we will fur- 
ther suppose that the field's velocity is such that it moves 
up the potential. If the change in the magnitude of the 
potential is negligible over a given cycle, this cycle can 
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FIG. 1: Illustrating the general features of the phase space 
discussed in the text that are relevant to a graceful entrance to 
inflation. The variables are {H, a} and all trajectories evolve 
in an ant i- clockwise direction. The top diagram corresponds 
to < V < Fd-it and the bottom panel to V > V C iit- In the 
former case, the important features for a graceful entrance 
are an area of phase space in which oscillations take place, as 
indicated by the region enclosed by the dotted line, and tra- 
jectories which evolve into an inflationary phase. The centre 
equilibrium point about which oscillations occur is marked 
with a circle, the separatrix is represented by the dotted line 
and the saddle equilibrium point occurs at the self-intersection 
of the separatrix. In the lower panel (V > V CI it), there ex- 
ist no centre or saddle equilibrium points and all trajectories 
evolve into an inflationary phase. 



still be represented as a trajectory in the phase space as- 
sociated with an instantaneous value of the potential, as 
illustrated in the top panel of Fig. ^ Over a large num- 
ber of cycles, however, this 'instantaneous phase space' 
will become significantly modified. Indeed, if the poten- 
tial energy increases monotonically as the field evolves, 
the region of phase space in which oscillations are pos- 
sible will become progressively smaller. Since the oscil- 
latory region disappears completely at V = V CT it, this 
corresponds to the largest value that the field's poten- 
tial energy can acquire before the oscillations are broken. 
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At this point the trajectory in the phase space rapidly 
evolves into the inflationary regime corresponding to the 
lower panel of Fig. ^ Hence, the value V cr it sets the 
energy scale at the onset of inflation. Once inflation be- 
gins, the field's kinetic energy will rapidly tend to zero 
as the field slows down, reaches a point of maximum dis- 
placement and rolls back down the potential. Slow-roll 
inflation will then arise if the potential satisfies the usual 
slow-roll constraints. 

Although the dynamics may be further complicated 
by the introduction of additional matter sources [19|, the 
above description outlines how a field can evolve up its 
self-interaction potential while the universe undergoes os- 
cillations. This forms the basis for a graceful entrance 
mechanism that naturally generates the conditions for 
slow-roll inflation. A phase plane analysis is important 
since it highlights the relevance of the centre and saddle 
equilibrium points. In particular, a sufficient but not nec- 
essary condition for cyclic behaviour is the existence of a 
centre equilibrium point, while a saddle point is required 
to separate those regions of phase space where cyclic and 
non-cyclic behaviour takes place. Furthermore, the dis- 
appearance of the saddle point at a critical value of the 
potential is key to the graceful entrance mechanism. 

In the following Section, therefore, we will consider the 
conditions for the existence of saddle and centre equilib- 
rium points in cosmologies described by a set of general- 
ized Friedmann equations. 



and Eqs. Q and |J2J fully determine the cosmic dynam- 
ics. Differentiating Eq. (|TJ with respect to cosmic time 
implies that 



37cff 



(4) 



and Eqs. ancl then describe a closed dynamical 
system, where the Friedmann equation Q represents a 
constraint. The equilibrium points of this system arise 
whenever d = a = 0, which implies that 



7eff(a cq ) = - , H(a cq ) = 



(5) 



The stability of the system ©-0 can be determined 
by linearising about the equilibrium points and evaluat- 
ing the corresponding eigenvalues. It is straightforward 
to show that the eigenvalues are given by 



A 2 = -; 



1 Cfyeff 

a 2 din a 



(6) 



In general, therefore, a necessary and sufficient condition 
for the equilibrium point to be a centre is 



A 2 < 



H = , 



7eff = g , 



dy, 



c rr 



din a 



>0 , 



whereas the corresponding condition to be a saddle is 



(7) 



III. GENERAL CONDITIONS FOR CENTRE 
AND SADDLE EQUILIBRIUM POINTS 

A. Relativistic Cosmology 

We wish to study classes of cosmological models that 
are motivated by the braneworld paradigm and to de- 
termine whether such models display the general char- 
acteristics required for a graceful entrance to inflation. 
However, we will first consider the relativistic cosmology 
based on classical Einstein gravity, since this system pro- 
vides a suitable framework for considering more general 
models. Specifically, we will consider a positively-curved 
FRW cosmology sourced by an effective perfect fluid with 
an energy density and pressure related by an arbitrary 
equation of state, p e g = [jeS — l]peff , where it is assumed 
implicitly that the equation of state parameter is a known 
function of the scale factor, y c s = 7eff(a)- 

The Friedmann and fluid equations for this model are 
given by 



H 2 



87^ 
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Peff 



PcS 



3 

= -3Hy cS p cS 



(1) 
(2) 



where 



A > 



H = 



7eff=3 



din a 



< 



(8) 

It is interesting to note that under the assumptions we 
have made, centre and saddle points are the only types 
of equilibrium points permitted. 



B. Braneworld Scenarios 

The four-dimensional cosmological dynamics for a wide 
class of positively-curved FRW braneworld scenarios can 
be modeled in terms of a generalized Friedmann equation 
of the form 



H 



/(«) 



(9) 



where p is the total energy density of the matter con- 
fined to the brane. The function L(p) is assumed to be 
positive-definite and parametrizes the departure of the 
model from the standard relativistic behaviour, H 2 oc p. 
The function /(a) is a function of the scale factor and, 
in a braneworld context, usually parametrizes the effects 
of a bulk black hole on the four-dimensional dynamics 
[32l l33l l34| . If the matter fields are confined to the brane, 
they satisfy the standard conservation equation 



d = Ha , 



(3) 



p = -3Hyp , 



(10) 
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where the equation of state parameter is defined by p — 
(7 — l)p and is once more to be viewed implicitly as a 
known function of the scale factor. 

The system l|9"))- H10|) can be expressed in the standard 
form of Eqs. Q~© by defining an effective energy den- 
sity m 



Peff = pL + 



/(«) 



(11) 



Differentiating Eq. Ijllfl and comparing with Eq. then 
implies that the effective equation of state parameter has 
the form 



7eff 



, r2 M 2d(L 2 ) 



1 



df 



87r£pj din a 



(12) 



In principle, therefore, if the conservation equation 
(110(1 can be solved for a given equation of state, 7(a), 
the effective equation of state (|12[1 will be a known func- 
tion of the scale factor once the functional forms of L(p) 
and f(a) have been specified for a particular braneworld 
model. Differentiation will then yield the necessary in- 
formation required in order to determine the existence 
(or not) of the centre and saddle equilibrium points, as 
summarized in Eqs. Q and Jjyl. 

In the following Section, we will employ these results 
within the context of a specific braneworld model. 



IV. SHTANOV-SAHNI BRANEWORLD 

The Shtanov-Sahni (S-S) braneworld scenario embeds 
a co-dimension one brane with a negative tension a in a 
five-dimensional Einstein space sourced by a positive cos- 
mological constant, where the fifth dimension is timclike 
[2fil l3fil ] . The effective Friedmann equation on the brane 
is given by [2(| |3|| 



H z = 



8w£ : 



pi 



,>2 1 



2(7 



771 
~7A 



1 



(13) 



where a = —a. If the bulk space is not conformally 
flat, the projection of the five-dimensional Weyl tensor 
induces an effective 'dark radiation' term on the brane, 
parametrized by the constant m [33 , l33L l3^ | . A confor- 
mally flat bulk corresponds to m — 0. Under quite gen- 
eral conditions, a negative quadratic dependence on the 
energy density in the Friedmann equation (|13fl implies 
that a collapsing brane world-volume is able to undergo 
a non-singular bounce |2r| . 

Comparison of the Friedmann equations 10 and i|13|) 
implies immediately that 



L 



1 



_P_ 

2a 



/(«) 



m 
~7a 



(14) 



and substituting Eq. I|14|) into Eq. I|12fl implies that the 
effective equation of state parameter is given by 



7eff 



7p(l - p/cr)a 4 + m/2irl 



Pi 



p(l - p/2a)a 4 + 3m/8ir£ [ 



(15) 



pi 



Since we are interested in whether a graceful entrance 
to inflation can occur in this model, we will determine the 
equilibrium points that arise when the matter confined 
to the brane corresponds to a scalar field that is rolling 
along a constant potential, V. This is equivalent to a 
matter source comprised of a massless scalar field and 
a cosmological constant. The equation of state for the 
system is given by 



7 



2 1 



(16) 



In general, the equilibrium points for this system will 
arise whenever Eqs. J7J or JSJl are satisfied. It follows 
from Eqs. I|13|) an d (|15fl that such points occur when 



2cr / a 4 a 1 
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Ant- 



(18) 
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For finite values of the scale factor, Eq. I|17|l may be 
simplified after substitution of Eq. 1|18|) : 



1 9 
— = A-Ki^p 



t_^P +1 (P_ l 

3 3(7 \<7 



(19) 



and Eq. I|19|) may then be employed to express Eq. 1|17|) 
in the form of a quartic equation in the energy density: 



5p; 



2A 



3(7 



V 



V 



G(7 



V 



Peq -v = o 



(20) 



The solutions to Eq. (|20|l yield the values of the energy 
density at the equilibrium points and the corresponding 
value of the scale factor can then be deduced directly 
from Eq. I|19|) . For physical solutions, one must ensure 
that the energy density and scale factor are positive and 
real at each equilibrium point. 

The dark radiation on the brane can significantly influ- 
ence the dynamics of the system. In view of this, we con- 
sider separately the cases where this radiation is present 
or absent in the following Sections. 



V. NO DARK RADIATION 

If no dark radiation is present (m = 0), the quartic 
equation l|20|) reduces to the quadratic constraint 



ioa 



3(7 



8 W 
3 + 



w = 



(21) 
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and this equation can be solved in terms of the field's 
kinetic energy: 



4ct - 4V ± 2 [W 2 - 18Va + 4cr 2 ] : 



(22) 



It follows that there can be at most two static (physi- 
cal) solutions to the Friedmann equations © and (|10fl . 
depending on the magnitude of the potential. If V < 0, 
there is only one solution to Eq. (|21|l with <\? > 0, im- 
plying there is only one equilibrium point in the phase 
space. For V = 0, there is one solution with </> 2 > 
and a second, but physically uninteresting, point where 
the field's kinetic energy vanishes and the scale factor di- 
verges. For a positive potential, on the other hand, there 
are two real roots to Eq. I|21|l if the condition 



W 1 - 18Va + 4a 2 > 



(23) 



is satisfied. The values of the potential which bound 
the region of parameter space in which there are no real 
solutions are therefore given by 



V = 




We find that for < V < V cr it , where 




(24) 



(25) 



there are two real solutions to Eq. J2J. Moreover the 
roots are physical, since they are positive and lead to real 
values of the scale factor. The roots merge and disap- 
pear at Krit- Although other solutions to Eq. (|21|l arise 
when the potential is greater than the positive branch of 
Eq. i|24[l. these do not correspond to physical equilibrium 
points, either because the solution requires <p 2 < or the 
value of the scale factor is imaginary. 

Hence, there are two static equilibrium points for < 
V < V CI it ■ The nature of these points can be determined 
from condition © after substituting Eqs. (|15fl . i|17|) and 
(|21|l and we find that the eigenvalues are given by 



A 2 = 47r4i<^ 



4- 



a V 



+ 4V 



(26) 



This implies that a given point will correspond to a centre 
if (f> 2 > | (<j — V), otherwise it will represent a saddle. 
Consequently, for physical roots to Eq. (|2~T|l . the field's 
kinetic energy has a value 



at the saddle. 

It follows from Eqs. Q27JI and Ij28(l that the centre and 
saddle equilibrium points move towards each other as the 
value of the potential is increased, eventually merging 
and disappearing at V — V CT i t - These are precisely the 
requirements for a graceful entrance to inflation that were 
outlined in Section II. 

We now proceed to illustrate the cosmic dynamics 
in the phase space. The dynamics is clearly two- 
dimensional due to the Friedmann constraint equation 
(113(1 . In principle, therefore, the phase space could be 
represented in two dimensions by introducing appropri- 
ate variables. However, this may be non-trivial in prac- 
tice, since more than one set of variables may be required. 
(This would then imply that more than one phase plane 
diagram would be needed in order to present the full 
dynamics). In the context of the present discussion, 
it proves convenient to numerically integrate the field 
equations and then view the evolution of the universe 
in the three-dimensional space spanned by {a,H, cf>}. 
Since the Friedmann constraint defines a surface in this 
phase space, the evolution of the universe can then be 
parametrized in terms of a trajectory on this surface. 
In the following, we will refer to such a surface as the 
Friedmann surface. Identifying a coordinate system that 
covers the Friedmann surface is then equivalent to finding 
variables in which the dynamical system becomes explic- 
itly two-dimensional, and the need for more than one 
coordinate patch to cover the surface is then equivalent 
to requiring more than one set of variables to complete 
the phase space. 

We focus on the case of a positive potential since we 
are interested in establishing the conditions for inflation 
via the graceful entrance mechanism. Fig. illustrates 
the Friedmann surface and the integrated trajectories on 
this surface when V < V CI it ■ The lower part of the figure 
represents the projection of these trajectories onto the 
two-dimensional plane spanned by the variables {(j), H}. 
The two equilibrium points representing the centre and 
saddle are shown. Fig. [3]illustrates the corresponding dy- 
namics for V > V cr it, where the equilibrium points have 
effectively merged. All trajectories now evolve through a 
bounce into an ever-expanding inflationary period of de 
Sitter expansion with H — ► constant, (f> — > and a — > oo. 

Comparison between the phase spaces of Fig. ^ and 
those of Figs. El and confirms that a graceful entrance 
to inflation, as outlined in section II, can occur in this 
model. 



VI. EFFECTS OF DARK RADIATION 



• 2 -4V + 4(7 + 2y/9V 2 - 18Va + 4ct 2 
^c q = jj » (27) 

at the centre equilibrium point, whereas it takes the value 



-AV + 4ct - 2y/9V 2 - 181/(7 + cr 2 



(28) 



In the case where dark radiation plays a dynamical 
role in the Friedmann equation (|13|) . the full quartic ex- 
pression given by Eq. (|2Tj|> must be solved to determine 
the nature of the equilibrium points. One must also en- 
sure that the scale factor and kinetic energy of the field 
take positive and real values at these points. The con- 
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FIG. 2: The top panel represents the phase space in the vari- 
ables {H, <j>, a}, when m — 0, a = 0.05 and < V < Vait, 
where V cr it is defined in Eq. I|25|l . Numerical values are given 
in Planck units. The surface defined by the Friedmann equa- 
tion with p — 4> 2 /2 + V, is also shown and all phase space 
trajectories lie on this surface. The compactified coordinates 
x = arctan(//), y = arctan(ln^) and z = arctan(lna) have 
been employed in order to show the entire phase space. The 
lower panel represents the two-dimensional projection of this 
space onto the plane spanned by the variables {H, <fi}. The 
centre equilibrium point is denoted by the solid circle and 
the separatrix is represented by the dotted line. The saddle 
point occurs at the point where the separatrix self-intersects. 
The axes have been compactified using the coordinate change 
x = arctan(//) and y = arctan(ln0). 




H 



FIG. 3: As for Fig. |5J but now for a potential V > Krit- 
There are no finite equilibrium points in the phase space in 
this region of parameter space. 



_ /V^ _ J_V = ^ _ JL (29) 
V a 2 2^ J 3 4m ' V ' 

Provided V < a, the hyperboloid's topology depends 
only the relative values of the dark radiation parameter, 
to, and the brane tension, a. Indeed, there is a criti- 
cal value of to at which the topology changes and this is 
given by 



straint 12UI) can be solved analytically using the standard 
techniques for quartic equations, but since the result- 
ing expressions for the roots are very lengthy, we will 
not present them here. However, further insight may be 
gained by investigating the topology of the Friedmann 
surface. This may be achieved by rewriting the Fried- 
mann constraint equation l|13|) in the form of a hyper- 
boloid: 




Consequently, the surface defined by Eq. I|29(l consists 
of two disconnected pieces if to < m cr it, whereas it is 
a single surface for to > TO cr it. It is possible that one 
of the two surfaces may correspond to a region where 
4> 2 < 0, in which case it would be unphysical. However, 
the qualitative dynamics of the universe will clearly be 
radically different above and below the critical value m cr i t 
and we therefore proceed to consider each regime in turn. 
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FIG. 4: Illustrating the position and nature of the equilib- 
rium points as a function of the field's kinetic and potential 
energies when < m < m cr it. A dashed line corresponds to 
a saddle point and a solid line to a centre point. The critical 
value V C rit2 represents the magnitude of the potential when 
one of the centre points and the saddle disappear from the 
phase space. The numerical values chosen for the plot were 
a = 0.05 and m = 1.0 and the axes are measured in Planck 
units. 



A. m < m cri t 

Since the analytic expressions for the solutions to the 
quartic equation (|20|) are not particularly illuminating, 
we have employed them graphically to illustrate how the 
nature of the equilibrium points depends on the field's 
kinetic and potential energies for specific values of the 
dark radiation parameter and the branc tension. Eq. 
(|19fl was also employed to verify that the solutions to Eq. 
(I20|l correspond to physical values of the scale factor, and 
their stability was determined from Eq. (JBJ. 

The results are shown in Fig. @|for m = 1 and a = 0.05, 
where the locations of the saddle and centre equilibrium 
points for given field energies are represented by the solid 
and dashed lines, respectively. The qualitative behaviour 
is the same for all m < m cr jt- There is a critical value 
of the potential, V a n2, such that there are three real 
roots to the quartic equation (|2"U|) when < V < V cl -i t 2- 
These correspond to three distinct equilibrium points, 
two of which are centres whereas the third is a saddle 
point. At V — V r cr it2, however, the saddle and one of 
the centre points merge and, for V > V cr it2, these points 
correspond to unphysical roots. Consequently, there is 
only one equilibrium point for V > V CT H2 and this is a 
centre. 

The qualitative evolution of the universe in the regimes 
V < V cl -n2 and V > V cr i t 2 is shown in Figs. El and El 
respectively, where the trajectories have been calculated 
by numerically integrating the field equations for the spe- 
cific choices m = 1 and a = 0.05. The topology of the 
Friedmann surface in the upper panels of these figures, 
together with Eq. Q29[l. implies that the entire phase 



space can be represented as a single two-dimensional plot 
by employing cylindrical polar coordinates. These are 
shown in the lower panels of the figures. 

The dynamics in the upper sectors of the Friedmann 
surface in Figs. EJandEJis qualitatively similar to the case 
where no dark radiation is present (m = 0). The lower 
sector of the Friedmann surface contains the new centre 
equilibrium point that arises when m ^ 0. This point is 
real for all values of the potential. Hence, as shown in 
Figs. El and El the universe remains trapped in an indef- 
inite cycle of expansion and contraction if it is initially 
located in the lower half of the Friedmann surface. This 
behaviour can be understood since a sufficiently large 
(positive) cosmological constant introduces a large neg- 
ative effective cosmological constant in the Friedmann 
equation (|13|l as a consequence of the quadratic term in 
the energy density. This will prevent the universe from 
entering a phase of accelerated inflationary expansion. 
Indeed, although is is not relevant for the graceful en- 
trance mechanism, it can be shown that if V exceeds yet 
another critical value, the upper region of the Friedmann 
surface is no longer physical since the field's kinetic en- 
ergy becomes negative. In effect, therefore, the size of 
the universe is bounded from above. 

By comparing the qualitative dynamics of the universe 
when no dark radiation is present (Figs. EJandEl) to that 
illustrated in the upper sectors of the Friedmann surfaces 
in Figs. El and El we may infer that a graceful entrance 
to inflation, as outlined in Section II, may in principle 
occur for suitable choices of initial conditions. 



B. m > m cri t 

We may adopt a similar approach for m > m cr ;t. Fig. 
^illustrates the position and nature of the physically rele- 
vant equilibrium points as a function of the field's kinetic 
and potential energies for the specific choices m = 1.4 
and a — 0.05. The qualitative nature of this plot re- 
mains unaltered for all m > merit- It follows from Fig. 
that there is only a single saddle equilibrium point when 
the magnitude of the potential is below a critical value 
Vcrit3 ■ There is then a finite range of values of the poten- 
tial, Vc r it3 < V < Vc r it4, for which there are no physical 
equilibrium points. For V > I4rit4j on the other hand, 
there exists both a centre and a saddle equilibrium point, 
and at still higher values of the potential, the saddle point 
disappears once more. However, the equilibrium points 
which occur above V^ r it4 are not relevant to the graceful 
entrance scenario. 

Figs. El and El illustrate the Friedmann surface and 
and the corresponding phase trajectories for V < V cl -it3 
and Vr; r it3 < V < V CT it4, respectively. In Fig. El the 
saddle point represents a divide between cyclic and infla- 
tionary behaviour, although it should be noted that the 
cyclic dynamics does not occur around a centre equilib- 
rium point in this case. For V > V cr it3, the saddle point 
disappears and all trajectories eventually evolve into an 
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FIG. 5: The top panel represents the phase space in the vari- 
ables {H,(f>, a}, when m = 1.0, <j = 0.05 and < V < Vcrit2. 
The surface defined by the Friedmann equation 113L with 
p = <p /2 + V, is also plotted. The phase space trajec- 
tories lie on this surface. The axes have been compacti- 
fied using the rescalings x = arctan(//), y — arctan(ln0) 
and z — arctan(lna). The lower panel illustrates a cor- 
responding two-dimensional phase space. It follows from 
Eq. 1291 that there exists an axis of symmetry parallel 
to the a-axis through the point H — 0, <f> 2 — 2a — 2V. 
Cylindrical polar coordinates can then be defined by us- 
ing this axis. Specifically, we define a — a, X = Rcos# 

and Y = Rsin^, where X = ^ - -^)), 

Y^,andR= - ' + f=f«= - . The 

two-dimensional plot is then presented in the {8, a} plane, 
where we have once more compactified the a-axis using y = 
arctan(lna). This diagram is essentially a projection of the 
top panel about the axis defined above, with the points at 
6 = 7r identified with those at 6 = —it. The shaded areas 
mark the unphysical regions in this coordinate system. The 
centre equilibrium points are identified by a solid circle and 
the separatrix is represented by a dotted line. The saddle 
point occurs at the point where the separatrix intersects with 
itself. 




FIG. 6: As in Fig. but now for a potential V > V C rit2- 
There are no saddle points and only one centre equilibrium 
point in the phase plane. 




V 



FIG. 7: As for Fig. 0] but now for the case where m > m cr it. 
Numerical values chosen for the parameters are m = 1.4 and 
a = 0.05 in Planck units. 



inflationary region, as shown in Fig. As was the case 
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FIG. 8: The top and bottom panels are as described for Fig. 
13 but now for m > m cr i t and < V < V" C rit3- Numerical 
values for the parameters are m — 1.4 and a = 0.05 in Planck 
units. Trajectories representing cyclic behaviour arise in this 
case but there are no centre equilibrium points and only a 
single saddle point is present. 



in the previous subsection, it can be shown that for a 
sufficiently large value of the potential, the Friedmann 
surface becomes modified in such a way that the sur- 
face is bounded to always lie below a critical value of 
the scale factor. This implies that the size of the uni- 
verse is bounded from above and such behaviour follows 
once more due to the presence of an effective negative 
cosmological constant in the Friedmann equation. 

The question of whether a graceful entrance to infla- 
tion is possible if m > TO cr it is more difficult to answer. 
The key features of a cyclic region and a saddle point in 
the phase space, which disappear at a critical value of the 
potential, are indeed present. However, there is a further 
complication. For m < m cr it, the cyclic dynamics is al- 
ways dominated by the field's kinetic energy, but this is 
no longer the case when to > m cr it. 

Further insight may be gained from Figs. |S]andE] The 
shaded gray areas (red in the online versions) represent 
the regions of phase space where the field's kinetic en- 



FIG. 9: The top and bottom panels are as described for Fig. 
HI but now for m > m cr it and Vcrit3 < V < V cr it4. Numerical 
values of the parameters are specified as m = 1.4 and a = 0.05 
in Planck units. There are no equilibrium points for this case. 



ergy is negative, with the boundaries corresponding to 
the limit cp 2 = 0. For a realistic (i.e. sufficiently flat but 
field-dependent) potential, this implies that on a trajec- 
tory which passes sufficiently close to these boundaries in 
the 'instantaneous phase space', the field's kinetic energy 
will be so small that the gradient of the potential will be- 
come significant. This will result in the kinetic energy of 
the field falling to zero and the field turning around on 
the potential before it has climbed sufficiently far up the 
potential to drive a successful phase of slow-roll inflation. 
In the previous cases, such behaviour was only possible in 
the region of parameter space relevant to inflation, or in 
a region disconnected from it, such as the lower section of 
the Friedmann surface. In the present case, however, the 
two regions of the previously disjointed Friedmann sur- 
face are connected. In effect, therefore, the turn around 
in the field may occur too soon for inflation to occur if 

TO > TO cri f 

More specifically, let us assume as before that we have 
a realistic inflationary potential, with a magnitude ini- 
tially in the range < V < V C rit3, and that the universe 
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begins in the oscillatory region of the phase space. As 
the universe oscillates, the field climbs its potential. In 
terms of Fig. |H1 the instantaneous phase space is there- 
fore altered, with the saddle point moving to successively 
smaller values of the scale factor, and the shaded (phys- 
ically forbidden) regions growing in size. Ultimately, a 
point will be attained when the universe's trajectory ei- 
ther moves outside the region of the oscillations or it 
passes too close to the boundary of the lower, shaded re- 
gion. If the former behaviour arises, the graceful entrance 
dynamics applies as in the previous scenarios discussed 
above. On the other hand, in the latter case, the field 
may turn around on the potential without inflation occur- 
ring, since the scale factor will still be small at this stage 
and consequently the curvature term in the Friedmann 
equation will be significant. Such behaviour is qualita- 
tively similar to the effect discussed in the LQC scenario 
when perfect fluid matter sources are introduced into the 
system The issue of whether a graceful entrance to 
inflation may occur is therefore sensitive to the initial 
conditions. 



VII. DISCUSSION 

In this paper, we have examined the criteria that a cos- 
mological model must satisfy in order for it to undergo 
a 'graceful entrance' into a phase of inflation, whereby 
a scalar field is able to move up its interaction potential 
whilst the universe undergoes a large number of oscilla- 
tory cycles. Our approach has been to consider an ideal- 
ized model, where the matter is comprised of a scalar field 
evolving along a constant potential. This system provides 
a good approximation to more realistic scenarios where 
the potential is field-dependent if the change in the po- 
tential is sufficiently small over a large number of cycles. 
It therefore provides us with a methodology for identify- 
ing whether a particular cosmological scenario will meet 
the necessary requirements for graceful entrance. 

The mechanism we have outlined is very generic. The 
important ingredients are that the phase space should ex- 
hibit both a centre and saddle equilibrium point when the 
magnitude of the potential V falls below a critical value 
U cr it- These points gradually move towards each other 
as V increases and eventually merge when V — V&it- 
Above this scale, the points should disappear from the 
phase space. From the physical viewpoint, this behaviour 
arises because the potential now dominates the field's ki- 
netic energy, V > (f> 2 , and consequently drives a phase of 
inflationary expansion. It follows that the critical value 
V[; r it sets the energy scale for the onset of inflation. 

We have developed a framework for investigating a gen- 
eral class of braneworld models in this context that are 
characterized by a Friedmann equation with an arbitrary 
dependence on the energy density. Such a class of mod- 



els can be represented as a standard, relativistic Fried- 
mann universe by reinterpreting the effective equation of 
state of the matter on the brane. As a concrete example, 
we focused on the Shtanov-Sahni braneworld p6|. where 
the extra bulk dimension is timelike. We found that the 
question of whether the braneworld oscillates or expands 
indefinitely is dependent on the field's kinetic and poten- 
tial energies, as well as the dark radiation and the brane 
tension, a. If the dark radiation is negligible, the phase 
space has a similar structure to that outlined above and 
a graceful entrance to inflation can therefore be realised, 
where the energy scale at the onset of inflation is given by 
U cr it ~ a. The dark radiation can have a significant effect 
on the dynamics, however, and the question of whether 
a graceful entrance is possible when this term is present 
is sensitive to the initial conditions. 

In conclusion, therefore, we have presented a non- 
singular, oscillating braneworld that can in principle ex- 
hibit a graceful entrance to inflation. We anticipate that 
such behaviour will apply for a wide range of collaps- 
ing braneworld scenarios that are able to undergo non- 
singular bounces (see, e.g., pTjp and, moreover, will not 
be strongly dependent on the precise form of the in- 
fiaton potential. Similar dynamics was recently shown 
to arise in models inspired by loop quantum cosmology 
[TiilrillTil and it is interesting that such behaviour is also 
possible in braneworld models motivated by string/M- 
theory. 

Finally, our analysis also provides the basis for realis- 
ing the recently proposed emergent universe scenario in 
a braneworld context [Isl Issl]. In this scenario, the ini- 
tial state of the universe is postulated to be either an 
Einstein static universe or one that is oscillating about 
such a solution in the infinite past. The potential has an 
asymptotically flat section, where V < U cr i t as <p — ► — oo, 
but subsequently increases above V cr - lt once the value of 
the field has increased beyond a certain value. If the field 
is initially located in the plateau region of the potential 
and moving in the appropriate direction, it will eventu- 
ally reach the section of the potential which rises above 
14 r it, thus initiating inflation. Within the framework of 
classical Einstein gravity, the emergent universe scenario 
suffers from severe fine-tuning, since the static universe 
corresponds to a saddle point in the phase space. How- 
ever, a mechanism which allows for a graceful entrance 
to inflation provides a natural realization of this scenario 
with less fine tuning and a greater degree of freedom in 
the choice of inflationary potential. 
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